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Internal relaxation time in immersed particulate materials. 



o 
o 

> 

o 



>> 
i 

q: 

o ■ 

C/3 ■ 
>> 

Ph. 



in 



On 
O 



P. Rognon and I. Eina^Q 

School of Civil Engineering, JOS, The University of Sydney, Sydney, New South Wales 2006, Australia. 

C. GajQ 

Matiere et Systemes Complexes, Universite Paris Diderot-Paris 7, CNRS, UMR 7057, 
Batiment Condorcet, Case courrier 7056, F-75205 PARIS cedex 13, France. 
(Dated: November 13, 2009) 

We study the dynamics of the solid to liquid transition for a model material made of elastic 
particles immersed in a viscous fluid. The interaction between particle surfaces includes their viscous 
lubrication, a sharp repulsion when they get closer than a tuned steric length and their elastic 
deflection induced by those two forces. We use Soft Dynamics to simulate the dynamics of this 
material when it experiences a step increase in the shear stress and a constant normal stress. We 
observe a long creep phase before a substantial flow eventually establishes. We find that the typical 
creep time relies on an internal relaxation process, namely the separation of two particles driven 
by the applied stress and resisted by the viscous friction. This mechanism should be relevant for 
granular pastes, living cells, emulsions and wet foams. 

PACS numbers: 82.70.-y, 83.80.Iz, 02.70.Ns, 



I. INTRODUCTION 

The rheological behavior of soft amorphous materials 
such as foams [l], 3] , emulsions 0, 0, HI , living cells @ and 
granular pastes 7], [f| is an active field of research with 
many practical issues, for these materials are broadly in- 
volved in geophysical events, industrial processes and bi- 
ological organisms. Besides, understanding their dynam- 
ics presents great theoretical challenges as their struc- 
ture, similar to that of glasses, is disordered, metastable 
and out of thermodynamical equilibrium [Io| . 

Perhaps the most salient property of these materials is 
their ability to deform like a solid or to flow like a liquid 
depending on the loading conditions. At the microscopic 
level, this solid-liquid transition relies on the particles 
ability to move with respect to each other, similar to 
the cage effect in molecular or colloidal systems [ill, E3 ■ 
However, large enough bubbles, droplets and grains do 
not show significant Brownian motion and their solid- 
liquid transition is thus not thermally activated. The 
driving force of the relative particle motion comes only 
from the external loading, while the resisting force comes 
from the particle interactions and/or deformations. This 
force balance defines internal dynamics with an associ- 
ated time-scale, which is the core of the thixotropy ef- 
fects in various materials [IH, [ID, 0, [HI, HI] and of 
the relaxation time evidenced by oscillatory shear exper- 
iments [13, EH ■ The interstitial liquid is pivotal to these 
internal dynamics, as it affects both the driving force and 
the resisting force that the particles experience. 

Immersed particulate materials support the stresses 
not only through their contact network (known as os- 
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motic or effective stresses) , but also by the continuous liq- 
uid (pore pressure). This notion of pore pressure is a cen- 
tral point of soil stability analysis (see [l9[ and reference 
therein) , and was shown to be responsible for a long d elay 
in the triggering of submarine granular avalanches [201 ] . 
The underlying mechanism involves the dilatancy needed 
to start shearing and the time required to fill the grow- 
ing pores with liquid. This filling time then relies on the 
degree of drainage of the shearing zone, in other word on 
how long the liquid takes to flow from one part of the ma- 
terial to another. Assuming that the granular skeleton 
does not deform, the flow of liquid through the pore net- 
work satisfies a Darcy law with some permeability [2p| . 

Besides the pore pressure, the liquid has a strong effect 
on the particle interactions as it lubricates close pairs of 
surfaces. The viscous friction between two smooth rigid 
spheres is a well known problem (see for instance [2l|). 
However, the presence of asperities on rough surfaces has 
been shown to strongly affect lubrication, and this topic 
is still an active field of research [2^, HH, Hj, [25[ . Between 
bubbles or droplets, the flow of liquid is even more com- 
plex, as the surfaces are soft and can stretch [26|, [2?| , they 
can be permeable p8j and exhibit some repulsive electro- 
static interactions leading to disjoncting pressure [29[ as 
well as adhesion 0, H3, HH • 

This paper presents a numerical investigation on the 
dynamics of the solid-liquid transition in a model im- 
mersed particulate material. Our motivation is to cap- 
ture the common characteristics of bubbles, droplets or 
solid grains: all of these particles can deform elastically 
and interact through lubricated contacts. Therefore, we 
do not address the flow of the liquid in the pore net- 
work or any of the others interactions mentioned above. 
Instead, here we focus on elastic spheres interacting via 
lubricated contacts with a steric repulsion when surfaces 
get very close. We simulate the dynamics of this material 
using Soft Dynamics, the discrete element method that 
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we introduced in Refs. [32j,l33|. We subject large samples 
to a constant normal (osmotic) stress and a (smoothed) 
stepwise increase of the shear stress. We then measure 
the delay before the onset of a steady material flow. Fi- 
nally, we show that this macroscopic delay reflects the 
time taken by two particles for separating, a motion 
driven by a tensile force that results from the applied 
stress and that is resisted by the normal viscous friction. 



II. INTERNAL DYNAMICS: INSIGHTS FROM 
DIMENSIONAL ANALYSIS 

Dimensional analysis is often used to define the typi- 
cal basic times that the system involves. Let us recall its 
outcome, as expressed in various studies. We are dealing 
with elastic particles of size d, Young's modulus E and 
mass to. They are immersed in a Newtonian fluid (vis- 
cosity 77) and confined by an osmotic normal stress a yy . 
Four time scales emerge from those quantities. 

The first two time scales correspond to the motion of 
a single particle which experiences a driving force a yy d 2 . 
For dry grains, this motion is resisted only by the grain 
inertia (mass to), while for particles with negligible mass 
immersed in a fluid of viscosity r], it is resisted merely 
by the Stokes drag. The respective corresponding times 
are the inertia! time U [34l . I3H [36| and the Stokes time 
tat Hill: 
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The other two time scales can be directly obtained by 
replacing the scale of stress o~ yy by the Young modulus 
E in each above time scale. This leads to the time t c 
for sound to propagate through a grain (which can also 
be linked to the binary collision time [39() and to the 
time tM for an elastically deformed particle to recover its 
shape when immersed in a viscous fluid, a time similar 
to that of a Maxwell relaxation |40, El 
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(3) 
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Tab. U gives an estimate of those time scales for elastic 
grains and for bubbles. For the particle deformation to 
remain moderate, the reduced normal stress k = -^ L 
must be small, and time scale t c (rcsp. t^j) is thus shorter 
than ti (resp. tst)- 

The above references clearly show the usefulness of 
such time scales, especially for expressing flow laws. 
However, none of them accounts for the details of the in- 
teraction of a particle with its neighbourhood. The aim 
of the present paper is to point out the importance of 
such a time scale when particles are clastic and interact 
via lubricated contacts. 



TABLE I: Order of magnitudes of the internal time scales 
in immersed particulate systems. Time scales from dimen- 
sional analysis (defined in section [XXJ and delay time scale A 
from lubrication interactions (computed through the present 
study) . Particles consist in grains or bubbles of size d = 1 mm 
immersed in water (viscosity 10 -3 Pas). The grain density 
is 2.10 3 kgm~ 3 and their Young modulus is 10 10 Pa. The 
density of bubbles is null and their effective Young modu- 
lus is of the order of E « 10 2 Pa. We consider an osmotic 
normal stress of 10 7 Pa for the grains and lCT 1 Pa for the 
bubbles in order to obtain the same reduced normal stress 
K = a yy 

IE = 10" 3 in both cases. The delay A is estimated 
according to Eq.© with a steric length of ho = 1CP 7 m, which 
corresponds to particles in the Hertz limit. 
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"The above delay A was based on Hertzian elasticity. Of course, 
bubbles and droplets do not display any bulk elasticity, yet they 
do recover their shape upon load reversal. For the present purpose, 
their effective Young modulus was estimated as the surface tension 
(1CT 1 Nrn" 1 ) divided by the bubble size: E ss 10 2 Pa (see [H for 
more details on droplet contact). 



III. INTERNAL DYNAMICS: PARTICLE 
INTERACTION 



We consider an ideal suspension made of elastic parti- 
cles in a Newtonian fluid. We focus on dense configura- 
tions, where the surface-to-surface distance h between a 
particle and its few closest neighbours is typically much 
smaller than the particle diameter d. Under such con- 
ditions, the particle-fluid-particle interaction is mostly 
governed by the liquid that is present within a small gap 
of thickness h and radius a with h -C a <C d (FigJTJt) . 
Accordingly, the force and torque transmitted from sur- 
face to surface within each pair of neighbouring particles 
are proportional to their relative velocity. For such a gap 
with a low aspect ratio (h <C a) in which the flow satisfies 
the conditions of the lubrication approximation, the pro- 
portionality coefficients (which we here call viscous fric- 
tions) can be expressed rather easily (see Appendix [5} 
as functions of the gap thickness h and size a. 

As a result of the viscous force and torque transmit- 
ted through the gap, the particle surfaces locally deform 
elastically. Hence, for each pair of close neighbours, the 
combination of the gap viscous response and the particle 
elasticities is analogous to a Maxwell element (see Fig|TJ 
part al) in which both the stiffness and the viscous fric- 
tion are non-linear (in the sense that they depend on 
the gap geometry) and tensorial (since they reflect the 
normal, tangential and twist modes, see Fig. [5]). 

The viscous friction term that represents the normal 
relative motion (Poiscuillc flow) turns out to be typically 
much higher than the corresponding terms for the tan- 
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FIG. 1: (Color online) (a) Particle-fluid-particle interaction: (al) the force F ela that deforms the surface elastically decomposes 
into the viscous force F v ' a and the direct particle-particle force F dlr . The size a of the transmitting force region is represented 
in two asymptotic limits: (a2) the Poiseuille limit for non-deflecting spheres, and (a3) the Hertz limit when the deflection 8 
is much larger than the gap width h. (b) Plane shear simulated using Soft Dynamics. 1000 particles are subjected to normal 
stress a yy and shear stress a xy in a periodic domain (Lees-Edwards boundary conditions fi^V The height of the cell fluctuates 
freely in order to keep a yy constant. 



gential sliding mode (Couette flow) and the twist mode. 
Thus, the normal relative motion should be responsible 
for the longest relaxation time. Let us therefore recall 
the expression of the normal viscous and elastic forces 
(see Appendix [A} : 



7 4 



3 a' 
-irrj—h 



F? la = aES 



(5) 
(0) 



where E is the Young modulus of the particle and 8 
the normal surface deflection. The size a of the region 
that transmits most of the force basically depends on 
the distance between surfaces and on the flattening in- 
duced by their elastic deflections. If the elastic deflec- 
tion is much smaller than the gap thickness, the surface 
shape is mostly spherical and the size a is given by the 
classical lubrication approximation between two spheres: 
a « \fdti (FigQ]a2). We refer to this asymptotic limit as 
the Poiseuille limit. Conversely, if the elastic deflection 
is much larger than the gap thickness, the particle sur- 
face elastically flattens over a region whose size is given 
by Hertz's theory: a ta ^/d6 (Fig[T]a3). We refer to this 
asymptotic limit as the Hertz limit. 

So far, two particles could not statically sustain any 
nonzero confining force, since the fluid would keep flow- 
ing out of the gap and the particle surfaces would cor- 
respondingly keep approaching although in a slower and 
slower motion. In practice, however, there is usually a 
typical distance hg below which particle surfaces interact 



not only via the fluid, but also through a direct particle- 
particle force. This interaction can have various physical 
origins, among which direct contacts between asperities 
of rough solid surfaces [H, [H, [24|, |2f| or repulsion be- 
tween charged or polymeric surfactant molecules at the 
bubble or droplet surface (29[. Our present focus is not 
necessarily to capture any of the details of such effects, 
but rather to represent their common feature: a strong 
steric repulsion that works on the particle surfaces in par- 
allel to the viscous force (FigQ]al). For that purpose, we 
chose the following power law: 



jpdir pconf 



(7) 



For a 1, this law captures our requirement for a 
strongly increasing repulsion as the gap becomes thinner 
than the steric length ho, and a quickly vanishing repul- 
sion for thicker gaps. Ideally, we would have liked to have 
an infinite a, but as this function becomes more strongly 
non-linear, the critical time step for accurate simulations 
increases tremendously. A value of a = 8 was found as a 
good compromise between these two requirements. The 
force scale parameter is chosen as F£ on f — a yy d 2 . In this 
way, the surface-to-surface distances stabilize around ho 
when the material is subjected to a confining stress o~ yy . 
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IV. COLLECTIVE RESPONSE TO A STEP 
SHEAR STRESS 

In this section, we use Soft Dynamics to study the re- 
sponse of the material to a sudden increase in shear stress 
in plane shear geometry, in the absence of gravity, with 
a fully periodic domain (Fig[Tj3). The principle of the 
method is detailed in Appendix [B] The typical experi- 
ment consists of two steps: (i) the preparation of a dense 
compressed sample and (ii) the shear by itself. 

First, particles are set in a loose and random config- 
uration. They are then subjected solely to a designated 
target normal stress o~ vv . This makes the material denser, 
with closer particles, while the height of the periodic cell 
decreases. The material eventually reaches a static state 
where the confining stress is supported by some pairs of 
particles whose surface to surface distances are on the 
order of the steric length ho ■ A correct simulation of this 
compression should include the long range multibody hy- 
drodynamics interaction, while the particles are still in 
a loose configuration. However, we do not focus on the 
compression process and use it only as a way to obtain 
a close-packed configuration. In the second step, which 
is the prime focus of our paper, we subject the configu- 
ration to a (smoothened) step in shear stress from zero 
to a final value a xy . We then measure the time evolu- 
tion of the nominal shear strain 7. More preciseley, the 
shear stress is continuously increased from zero to its fi- 
nal value within a short but finite time of the order of tst 
(EqUJ) as a xy (l-e~~i^) (see Fig. El . Most of the results 
are obtained with (3=1. Different values ((3 = 0.5 and 
j3 = 2) will be used to probe the influence of the rate of 
the initial shear stress increase. 

We perform a scries of such numerical experiments for 
materials with various steric lengths ho, subjected to var- 
ious levels of shear stress o~ xy . The parameters of the 
system are summarized in Tab. [ill expressed in units of 
particle size d for the length, normal stress a yy for the 
stress and Stokes time tst for the time. The grains have 
a mass, but we set the inertial time much smaller than 
the Stokes time, so that the dynamics of the system is 
mostly driven by viscous interactions. The reduced nor- 
mal stress k = <j yy /E applied to the particles is much 
smaller than unity so that the elastic deformations would 
be much smaller than the particle size. We complete the 
current section by describing the resulting collective be- 
havior in our Soft Dynamics simulations, which we will 
then analyse in terms of the internal dynamics in Sec- 
tion El 



TABLE II: Properties of the simulated systems: inertial time 
U, reduced applied normal stress k, steric distance ho and 
shear stress a xy . 
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material first deforms very slowly. After this creep phase, 
it starts to flow steadily. During both the creep and the 
flow phases, j(t) is mostly linear (Fig. Eb)- Except for 
the very begining of the creeping phase, this behavior is 
not sensitive to the typical time of increase of the stress 
(Fig. El 3 , inset). We define the typical duration A of 
the creep phase using a geometrical construction that 
searches for the intersection between these two lines, as 
plotted in linear coordinates. Note that prior to the creep 
phase, there is another very short phase, during which the 
shear deformation quickly increases up to a small level. 
This corresponds to the immediate clastic deformation of 
the particles induced by the increase in shear stress. 

Fig. [3] shows the value of A for materials with differ- 
ent steric lengths ho- The creep duration A naturally 
increases when ho gets smaller. However, as long as the 
steric length is larger than roughly 10~ 2 d, A is small: 
only few times longer than the Stokes time. In contrast, 
it increases strongly for smaller ho- 



B. Varying the step shear stress 

In the second test, a material with a steric length ho = 
2.10 _1 d was subjected to varying step shear stress values 
(Fig. [H. It appears that below a critical shear stress 
threshold, the material creeps and then does not seem 
to flow, but rather reaches a static equilibrium within a 
finite shear strain smaller than 100%. Above this shear 
stress threshold, here found to be roughly 
the material creeps and then flows. One might imagine 
that the response delay in the creep phase before the 
onset of flow should be shorter if the applied constant 
shear stress is higher. Surprisingly, as can be seen from 
curves D-G of Fig. 01 i.e., with applied shear stresses 
that differ by a factor of 5, it seems that the duration of 
the creep phase remains essentially constant although the 
transition from creep to flow becomes sharper for larger 
values of the applied shear stress. 



A. The role of the steric length 

In the first test we subject materials with differ- 
ent steric lengths ho to a similar constant shear stress 
o~xij = 0.5o- yy (Fig. Ek)- For the longest steric distance, 
ho/d = 10 _1 , the shear strain increases almost steadily. 
In contrast, for the shortest distance ho/d = 10~ 3 , the 



V. MICROSCOPIC ORIGIN OF THE 
OBSERVED RELAXATION TIME 

Let us now interpret the macroscopic behaviors ob- 
served in the previous section in terms of the dynamics 
of two particles in relative motion. When the material 
is sheared, neighbouring particles need to move towards 
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FIG. 2: (Color online) Time evolution of the shear strain 7 induced by a rapid increase in shear stress a xy from to 0.5a yy 
between t = and i = tst- (a) Materials with a different value of the steric length ho. (b) Procedure to determine the typical 
duration A of the creeping phase, here for ho/d — 2.iCP' ! ; the dashed lines correspond to the linear fit of y(t) in the creep 
phase and in the flow phase; (Inset) for the same matrial, different rates of initial stress increase: j3 — 0.1 (black), 0.5 (red) 
and 2 (blue) (see introduction of Section H*VT) . 
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FIG. 3: (Color online) Macroscopic relaxation time A, esti- 
mated as shown on Fig. &jp, as a function of the steric length 
ho- Circles denote results from simulations, while the curve 
corresponds to 3 times the separation time t sep for two grains 
as denned by Eq. (S). The third data point (ho/d = 2.I(T 3 ) 
corresponds to curve F in Fig. 



and away from each other. This relative motion includes 
sliding, rotation and normal approach and separation. As 
shown in Appendix [XI the normal viscous friction is the 
highest, which implies that this mode of motion should be 
the slowest. As a result, the material lies most of the time 
in a situation where interparticle contacts can be viewed 
as frictionless as far as sliding is concerned, with only 
normal forces acting to bring particles into contact or 
away from each other in a slow motion. Hence, we focus 



FIG. 4: (Color online) Time evolution of the shear strain 7 
induced by establishing (from t = to t = tst) a constant 
shear stress (for larger i). Each curve corresponds to a dif- 



ferent value of the final shear stress: a xy /a y 



0.05 (A), 



0.1 (B), 0.15 (C), 0.2 (D), 0.3 (E), 0.5 (F) and 1.0 (G). Red 
(resp. black) curves denote tests in which the material even- 
tually flows (resp. does not flow). The same material (with 
ho = 2.10~ 3 ) was used for each run. Curve F corresponds to 
one of the experiments considered in Fig. [3] 



on the sole normal mode for interpreting the macroscopic 
observations. 

The duration of the normal separation process depends 
on the initial distance between surfaces, on the viscous 
friction and on the typical force transmitted by the fluid 
in the gap. The typical initial distance between the sur- 
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faces can be estimated as the steric distance ho. As the 
step in shear stress is applied to the dense material, the 
initial viscous friction depends on whether the interac- 
tion is initially in the Poiseuille regime (5 < ho) or in the 
Hertz limit (<5 > ho). Initially, the typical force transmit- 
ted within a pair of particles scales with a yy d 2 . Hence, 
the interaction lies in the Hertz limit when h /d < ni . 



Let us assume that, when the shear stress is applied, 
the typical transmitted force is still of the order of <J yy d 2 . 
Then, the time t sep for two particles to separate is such 
that J ( *_ a yy d 2 dt = f h l e £ o ^Ttrj^jdh, where the dimen- 
sion of the force-transmitting region is a = \/dh in the 
Poiseuille limit and a ~ \J d\8\ in the Hertz limit. Hence: 



in 



(S3 



£1 - 1 



±ln 



dK 3 



(Poiseuille, !f > 
(Hertz, ^«k3). 



(8) 



An arbitrary value of h sep has to be defined at the end 
of the separation process. Basically, it should be of the 
order of d. But the time of separation is only slightly 
depending on it. In the Poiseuille limit, the time of sep- 
aration is at the order of the stokes time, with almost no 
effect by the small steric length ho- In the Hertz limit, 
in contrast, it is strongly dependent on ho as well as on 
the reduced normall stress k = . When applied to the 
system simulated in the previous section (k = 10 -3 ), this 
simple separation time t sep captures the scaling behavior 
of the creep time A (FigJ3]) . The data are quantitatively 
well represented by the function bt sep involving a fit pa- 
rameter b ss 3. This constant is likely to rely on the de- 
tails of the contact network (number of coordination, het- 
erogeneities, anisotropy. . . ) and may also slightly depend 
on the geometrical definition of A (FigfSJa). The salient 
result is that, surprisingly enough, this very simple model 
with only one interaction and a single force scale o~ yy d 2 
captures the macroscopic delay observed: both t sep and 
A scales with ln/ip 1 while the interaction is initially in 
the Poiseuille regime (h > ft 3 = 10~ 2 ), and with h^ 2 if 
the interaction is initially in the Hertz regime (sec Fig. [3]). 

However, this model does not account for the sec- 
ondary effect of the shear stress we observed in Fig. [4] 
Basically, one could expect that under a larger ap- 
plied shear stress, the typical tensile pairwise normal 
force would increase. Accordingly, one should observe a 
shorter separation time. Fig. [4] reveals that the creep- 
ing time appears to be essentially independent of the 
applied shear stress. This may have some connection, 
although probably not direct, with the counterintuitive 
effect highlighted in (3^1 concerning the crossover be- 
tween the Poiseuille and the Hertz regimes: there exists 
an optimal tensile force that minimises the time of sep- 
aration between two particles. This optimal force cor- 
responds to the transition between the Hertz and the 
Poiseuille regimes. Thus, it depends on the gap h as: 

F opt w Ed 2 (see UK for more details). Increasing 

the force then brings the system further into the Hertz 
regime, which strongly increases the viscous friction and 



eventually results in a slower separation. 



VI. CONCLUSION 

The distinction between the flowing and non-flowing 
states of many materials is critical for determining fail- 
ure. The current paper focuses on studying the microme- 
chanical origins of this phenomenon in relation to ideal 
materials of closed-packed elastic particles in a viscous 
fluid. Using Soft Dynamics simulations, we show that 
when subjected to a step shear stress, this material can 
exhibit a long creep phase before a steady flow state sets 
in. The typical creep time seems to rely on an internal 
relaxation process: the normal separation of two parti- 
cles, which is driven by the applied stress and resisted by 
the viscous friction. This internal relaxation time thus 
depends on (i) the applied stress, (ii) the particle stiff- 
ness and (Hi) the steric length ho below which particle 
surfaces interact directly. 

The direct surface-to-surface interaction enables the 
material to statically sustain a confining stress, and 
whether or not the steric length is shorter than the typi- 
cal elastic deflection strongly affects the scaling behavior 
of the viscous friction. On one hand, the relaxation time 
is not very long in the Poiseuille regime (large ho and/or 
very stiff particles) as its scales with ln(/ig 1 ) and is thus 
only slightly longer than the Stokes time. On the other 
hand, the relaxation time can be very long in the Hertz 
limit (small ho and/or soft particles) as it then scales 
with ho~ 2 (see the estimation in table [XJ> . 

The long aging process corresponding to the Hertz 
limit should be a common characteristic of various im- 
mersed particulate materials like granular pastes, emul- 
sions and foams. We expect that the soft bubbles and 
droplets would reach the Hertz regime for relatively low 
confining stresses. For the solid grains, a rough estimate 
can be performed according to Hertz's elasticity which 

2 

defines the Hertz limit for ho < d (^jff-) 3 ■ Assuming that 
the steric length ho is determined by a surface rough- 
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ness of amplitude 0.1 /xm, that the grain size is d ~ 1 mm 
and that their Young modulus is E ~ 10 10 Pa, then the 
normal stress required to reach the Hertz regime would 
be on the order of a yy ~ 10 4 Pa. Note, however, that 
these materials involve further modes of interaction such 
as solid friction, interfacial rheology and adhesion, which 
would (i) induce other time scales and (ii) affect the hy- 
drodynamics within the gap and thus the expression of 
the viscous friction which the aging process relies on. In 
future works, we hope to include some such detailed in- 
teractions within our methodology. 
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APPENDIX A: LUBRICATION BETWEEN 
ELASTIC SURFACES 

The viscous flow between between two close elastic 
spheres is not a simple problem. However, simple scaling 
formulae can be deduced considering that only a small 
portion of the surface transmits the force (size a), and 
that its radius of curvature is large (a -C R). Then, 
the hydrodynamical interaction is similar to that which 
develops between two flat parallel disks of radius a sep- 
arated by h. The low aspect ratio of the gap, - C 1 
favors flows in the Stokes regime. 

A simple way to get the scaling formulae for the viscous 
force and torque is to consider the average shear strain 
rate e in the liquid, which scales like %- for a Couette 



flow and like 



3 V s 7ra 



for a Poiscuillc flow. On the one 



2 h 2iza h 

hand, the power dissipated by the viscous flow in a gap of 
a volume ira 2 h scales with ne 2 a 2 h. On the other hand, 
this power is equal to F ms v s for a relative velocity v s 
and to T ms w s for a relative angular velocity w s . This 
leads to the expressions of the viscous forces and torques 
in Fig. [5] Note that the viscous force and torque can be 
expressed through the viscous friction tensors Z and w: 



F ms = Z -V s 
~T ms = W ■ w 1 



with Z = £ n <X> n + X(I — n ® n) 
with W = a 2 Z, 

(Al) 

where n is the unit vector normal to the surface, / the 
unity tensor and £g> the standard cross product. 

The elastic forces and torques can also be estimated in 
a simple way by following the Hertz theory: the elastic 
strain e concerns mostly a small region of volume a 3 . It 
scales according to - for both the normal and the tan- 
gential deflections represented by S, and as — for a twist 



angle 9. The total elastic energy, Ee 2 , therefore scales 
like Ea 3 (f ) for deflect ions and according to Ea?9 2 for 
twists. The differentiation of these energies with respect 
to 5 and 9 gives the scaling for the elastic forces and 
torques expressed in Fig. [5l where all numerical constants 
are set to unity for the sake of simplicity. 



APPENDIX B: SOFT DYNAMICS METHOD 

We have developed the Soft Dynamics method to simu- 
late large samples of closel y p acked particles with elasto- 
lubricated interaction [32l. l33j| . In the present work, the 
particles are inertial so that their motion is deduced from 
the forces and torques they experience according to New- 
ton's law. as in the standard discrete clement methods: 



m i X l 

J i i 



. j? ext 
~p ext 



{ r t 



* x F el 

y x r ij 



(Bl) 



mi and Ji are respectively the mass and the inertia mo- 
ment of particle i, Xi and f2 j its position and its rota- 
tion velocity. Ff xt and T? xt are possible external force 
and torque that the particle experiences, which are both 
taken as zero in the present study. The sums run on each 
particle j interacting with i. The force Ffj a and the 
torque V fj a come from the elastic deflection and twist 
of the particle surfaces, and r y is the vector from the 
center of i to the interaction point. 

The uniqueness of the Soft Dynamics approach is to 
include the dynamics of the elastic deflection S ij and 
twist angle 9 ^ so that the force and torque balances are 



satisfied for each interaction. 



8 ^ and 



can be de- 



duced from the difference between the particle velocities, 

Xi and f2 ;, and the relative translational and rotational 
velocities of particle surfaces, v f • and wfj. Compatibil- 
ity then requires to have: 



ij = X j +n J x r 3 



Xi + Hi x r ■ 



1. 1 



(B2) 



For each of the interactions, the relative velocities v : 
and iuf • are given by the force/torque balances: 



^ dir 



p ela vis \ j^c 

ij ~ f\ij iQ 

p ela p vis _i_ p dir 

ij — ij ' ij ' 



(B3) 



and the fact that the viscous force/torque are propor- 
tional to them through the friction tensors Zij and Wij 
(see App. [£}: 



8 




F da = aE5 F da = aES T ela = a 3 E6 T ela = a 3 E6 



FIG. 5: (Color online) Scheme of the elasto-hydrodynamical interaction implemented in Soft Dynamics, including four modes 
of relative motion between particles with a surface-to-surface distance of h. Distance a is the size of the transmitting force 
region (see FigfjJ,). The blue arrows designate the flow velocity profile and the red spring the elastic deformation of the bulk 
material. The corresponding expression of viscous and elastic forces/torques are as follows, (a) The normal approach induces 
a Poiseuille flow and an associated viscous friction £ (also valid for a normal separation), (b) The tangential sliding induces a 
Couette flow and an associated viscous friction A. (c) Normal rotation (locally, the flow is like a Poiseuille flow), (d) Parallel 
rotation (locally, the flow is like a Couette flow). 



ij- a ■ I j j J (B4) 

ij ~ ij ' [ ij ij 

The second order equations (|B2[) are integrated using 
a standard predictor-corrector scheme, while the first or- 
der equations (|B3[) are integrated using a simple Euler 
scheme. The time step of integration is a fraction of 
the shortest time scale in the system, namely tm = 
Through careful examination we found that a time step 
of five time smaller than tM is short enough. 



Since we do not include long range hydrodynamical in- 
teractions, it is sensible to introduce a cut-off distance: 
here two particle interact only if their surface to surface 
distance is smaller than 0.25d. The size of the trans- 
mitting force region a is then being interpolated between 
the two asymptotic limits of the Poiseuille and the Hertz 
models, as a way to get a continuous transition. The 
interpolation we choose is: a = ^d{h+ \8\). Note that 
different choices of interpolation are possible, although 
this should have no incidence except in the vicinity of 
the crossover between both scaling regimes. 
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